Point defect in solids: Shear dominance of the far-field energy 
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C^ ■ 

^-^ , It is shown that the elastic energy far from a point defect in an isotropic solid is mainly shear 

^— «, ■ elastic energy. The calculation, which is based on a standard dipole expansion, shows that no matter 

^vj ' how large or small the bulk modulus is compared to the shear modulus, less than 10% of the distant 

point defect energy is associated with volume changes. 
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PACS numbers: 61.72.Ji 



An isotropic solid has two independent elastic constants, the shear modulus G and the bulk modulus K. It has 

been suggested that when a physical property depends on both G and K, the dependence on the shear modulus is 

^ often the most important ^]. Examples of this "shear dominance" were given by Granato in his important paper 

from 1992 [2]; they include the fact that defect energies vary only slightly with the bulk modulus, but arc linearly 

dependent on the shear modulus (for Cu he estimated that only about 3% of the energy of an interstitial is bulk elastic 

energy). In Ref. [l|, which dealt with the non-Arrhenius viscosity of viscous liquids and its possible explanation in 

terms of temperature-dependent instantaneous elastic constants [3], the following result was proved: If temperature 

dependence is quantified in terms of log-log derivatives, at least 92% of the temperature dependence of the molecular 

vibrational mean-square displacement over temperature comes from the instantaneous shear modulus, whereas at 

'T' ' most 8% is due to that of the instantaneous bulk modulus. 

'^ , How general is "shear dominance"? Is it coincidental, or is it a fairly general theme of three-dimensional elasticity? 

G Of course, a property like the Poisson ratio depends equally on the bulk and shear moduli, so shear dominance cannot 

be a completely general phenomenon. Nevertheless, there are several examples of it so it deserves to be investigated. 

As one contribution to shedding light on this question, we below calculate the maximum ratio of bulk elastic energy to 

shear elastic energy far from an arbitrary point defect in an isotropic solid. Based on a standard dipole expansion [J| 

it is shown that less than 10% of the distant elastic energy is bulk elastic energy, i.e., associated with density changes. 

Although our main objective is to identify the dominant contribution to the elastic energy far from a point defect, 

it should be noted that according to Zener's "strain energy model" [5| most energy associated with a point defect 

t^^ , is elastic. Elastic energies of point defects fall of rapidly (as r~^), as evident from Eshelby's famous solution of the 

^vO ■ problem of an elipsoidal inclusion '(i lD, [sl, l9|, llO|, [ill, [l^ , a result which is general. This means, of course, that most 

of the elastic energy is located close to the point defect and that a dipole expansion is not realistic for calculating 

the total defect energy. Nevertheless, our results may be taken as an indication of what contributes most to the total 

defect energy, the shear or the bulk elastic energies in the defect surroundings. 

C^ ' The nature of the far-field deformation is also important for understanding and modeling the long-range properties 

of defect-defect interactions. If, for instance, the displacement field is dominated by shear displacements, a defect will 

I ■ only interact weakly with one far away that is dilational close to its center. Some time ago Andreev discussed the 

' '^ unique topological characteristics of point defects in three dimensions 13 1. A comprehensive review of point defect 

b^ ' properties focusing on the vacancy was given by Kraftmakher 14 1 . The question of the nature of the deformation 

Q [ far from a point defect is relevant also for applications in materials science, e.g., for understanding fracture [9|]. In 

1!V ' its macroscopic description the question of the far-field properties of a point defect has been studied because of its 

. ^H [ important for understanding the mechanical properties of polycrystals and composites [J, lo, ll5| . In this contex Onaka 

^ ' recently calculated the elastic shear energy due to a macroscopic point defect (an "inclusion") |l6|, but did not study 

M ] the general case, leaving undetermined what is our focus here, the ratio between shear and bulk energy (see also Ref. 



[12|). Garikipati and co-workers [17[ recently discussed the role of continuum elasticity in determining the formation 
energy of a point defect utilizing Eshelby's result for the work done by an external stress during the transformation 
of an inclusion [6| , an interesting paper that also elicidates the limitations of this approach as compared to atomistic 
calculations. We mention these works also to emphasize the continuing interest in the far- field properties of point 
defects and their macroscopic analogs. 

First some preliminaries: We use the standard Einstein summation convention. The point defect is located at the 
origin, is modelled as follows: Imagine a small sphere of radius R surrounding the point defect with all atoms within 
the sphere removed. The effect of the defect on the surroundings are represented by a suitable distribution of external 
forces Fi applied to the surface of the sphere. We define a tensor Ay as the following integral over the surface of the 



sphere (where dA is the area element): 

Ay = / F,iy)yjdA. (1) 

J\y\=R 

Because the force distribution models the effect of the atoms within the sphere on the surroundings, the angular 
momentum of the force distribution must be zero. This implies that A is symmetric: Xij = Xji. 

Sharp brackets {)r, referred to as "averages," denote integrations over the surface of a sphere with radius r ^^ R 
centred at the origin. The following identities become useful later on: 



{XiXjiT- — ^ Oij yZij 



and 



{xiXjXkXi)r = — [SijSu + SikSji + SiiSjkj ■ (3) 

If (Xxx) is an abbreviation of XijXiXj, we find from these identities that 



{iXxx))r = ytrA (4) 



and 
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Moreover, if Ax is the vector whose i'th component is XijXj, we have 



{{Xxxf)r^'^({trXf + 2tv{X^)). (5) 



((A^xx)}, = ((Ax)2}, = -tr(A2). (6) 

As mentioned, far from the defect both the bulk and the shear elastic energies vary with the distance from the 
origin as r~^ . When one averages over a sphere with radius r ^ R, the result for both the bulk and the shear elastic 
energy must be a scalar varying with distance as r~^ that is of second order in the forces Fi. It follows from the below 
calculation that these two scalars are both uniquely determined by G, K, and the A-matrix. Consequently, because 
these two scalar functions are of second order in the forces, the ratio of bulk to shear elastic energy must have the 
following general structure: 

Bulk elastic energy ^(trA)^ + i3tr(A^) 

Shear elastic energy ^ C(trA)2 + Dtr{X'^) ' ^ ' 

Defining a — tr(A^)/(trA)^ this ratio is {A+ Ba) / {C + Da) . The quantity a varies between 1/3 and oo: By normalizing 
we may assume that trA = 1; if the eigenvalues are denoted by /i^ we thus have /ii + /i2 + Ms = 1 and consequently 
a = fj-l + fj,2 + ("i- From this it follows that a varies between 1/3 and oo. Since the energy ratio of Eq. ([7]) is 
a monotonous function of a and since perfect isotropy (a — 1/3) implies zero bulk elastic energy (in this case the 
displacement field is radially symmetric varying as r~^ which implies no volume changes), the maximum bulk to shear 
elastic energy ratio arises in the limit a — > oo. Thus in the below calculation we may ignore all terms with trA and 
keep only terms with a tr(A^) factor. 
Poisson's ratio a is defined [18| by 

1 3K-2G ,^, 



If a force F is applied at the origin of an isotropic elastic continuum, the displacement field at the point {x' , y' , z') is 



given [18j by 



, , , M (3 - 4a)F + (F ■ n>^ 
u{x ,y ,z) ex , (9) 

where r'^ — x'"^ + y'"^ + 2'^ and n' — {x' ,y\z')/r' is the unit vector pointing from the origin to {x',y',z'). It is 
convenient to introduce the variable 

A = 2-4cr, (10) 

in terms of which Eq. ([9]) becomes 

u, oc (A + l)^ + ^-p. (11) 

We proceed to perform a standard dipole expansion by first noting that, if y is the coordinate for a point on the 
small sphere surrounding the defect (|y[ = R) and x is the coordinate for the point of interest far away, to lowest 
order in |y|/|x| — i?/|x| we have if r = |x| 



|x-y| 

Similarly 



|x-yP 



x» + y^-2x.y)-'^ = .-(l-22^)-" ^.-H-i:^. ,12) 



/ 2 , 2 o N-3/2 -3/1 oX-y\-3/2 x-y 

= (x +y -2x-y) =r (1-2^^ =r +3 — =- . (13) 



To calculate the displacement field at point x we first note that when the force -Fi(y) is integrated over the small 
sphere radius R, the result is zero. Thus when one integrates over the small sphere, the first term of Eq. pip to 
lowest order in l/r becomes 

(A + 1) / lM^dA = {A + l)i F{y)^dA 
J|y|=fl|x-y| J|y|=R, r3 

= (A + 1) r-^KjXj . (14) 

Similarly, to lowest order the second term of Eq. pip gives the following contribution to the displacement field 

-^ '- '- dA ^ f Fj{y)[r •" + 3—r-){xj -yj){xi~y^)dA 

r^ XjiXj — r^ XjjXi + 3r^ XjiXjXiXi . (15) 



\y\=R |x-y|3 J|y|=;j^ "V r 



Thus, if proportionality is replaced by equality for simplicity of notation - which is OK because we only wish to 
calculate an energy ratio and have already dropped the overall proportionality constant of Eq. ([H]) - we find 

Ui — Kr XijXj — r~ trA xi + ir XjiXjXiXi . (16) 

Next we calculate the stain tensor. First, note that if dk is the partial derivative with respect to Xk, we have 
dkr-"^ = (-n)r-("+2)a;fe. Thus 

dkUi = —3Ar~^XijXjXk + Ar^^XijSjk + 3r~^tTXxiXk — r^^trXSik — Ibr^"^ XjiXjXiXiXk 

+ 3r^^XjiSjkXiXi + 3r~^XjiXjSikXt + 3r~^^XjiXjXi6ik . (17) 



All terms in this expression vary with r as r~^. Consequently, this term is common to the bulk and shear elastic 
energies and may be dropped from our calculation of their ratio. The calculation is simplified notationally by putting 
r — I (or, equivalently, replacing Xi by Xi/r). When this convention is adopted, the strain tensor is given by (ignoring 
the factor 2 in the strain tensor definition, Uik = {diUf. + dkUi)/2) 

Uik = 2 (— trA + 3{Xxx)) dik + 2AXik + GtrXxiXk — 30{Xxx)xiXk 

= 2 (—trA + 3{Xxx)) dik + 2AAifc + (6trA — 30(Axa:)) XiXk 

+ (6 - 3A) {XijXjXk + XkjXjXi) (18) 

Throwing out terms with trA we end with the following expression for the strain tensor 

Uik = 6{Xxx)Sik + 2AXik - 30{Xxx)xtXk + (6 - 3A) (X^jXjXk + XkjXjX,) . (19) 

To calculate the bulk elastic energy we need the trace of this which, when again terms with trA are dropped, is 
given as (xiXi = r^ is put equal to unity) 

trw = 18(Axa;) - 30(Aa;a;) + (6 - 3A)2(Aa;a;) = -6A(Aa;a;) . (20) 

The bulk elastic energy density [18| averaged over the sphere with radius r (subsequently put equal to unity) is, when 
use is made of Eq. ([5]) and terms with trA are ignored, given by 

Bulk elastic energy = — {{tT{u))^)r = — 36A2((Aa;x)2),, = —KAhriX^) . (21) 

2 2 5 

To find the shear elastic energy we need the transverse part of the strain tensor, w^. (in terms of which the shear 
energy density is Gw^^i) defined as [l^: [1^ 

■"ifc = Uik - ■::tr(u)Sik 

= 2(3 + A){Xxx)Sik + 2AXik — 30{Xxx)xiXk + 3(2 — A) {XijXjXk + XkjXjXi) . (22) 

Squaring and summing of all the elements of the transverse strain tensor, which is required to calculate the shear 
elastic energy, leads to the following (xiXi = r^ is again put equal to unity): 

uji^ui = 4(3 + AfiXxxfS + AAhr{X^) + 900(Aa;a;)2 + 9(2 - A)2[2(Ax)2 + 2{Xxxf] 

+ 8A(3 + A)(Axx)trA - 120(3 + A)(Axa;)^ + 12(3 + A)(2 - A)2{Xxx)^ - 120A(Axa;)^ 

+ 12A(2-A)2(A2xa;) - 180(2- A)(Aa;2:)2(Aa;a;). (23) 

Averaging this expression (ignoring all trA-terms) leads to: 

{uikuik)r = (l2(3 + A)2A+4A2 + 900^ + 18(2-A)2(i + A) 

- 120(3 + A)— + 24(3 + A)(2- A) 120A— 

15 15 15 

+ 24A(2-A)i-360(2-A)— ')tr(A2) 
3 15/ 

= - (7A2 + 12A+108)tr(A2). (24) 

Summarizing, we find that 



Bulk elastic energy (12/5)KA2 _ K QK^ 

— /o l^\r^(n A a I 1 o A I TTvoT ^ 7-<' -7 A 2 I 1 o A I 1 no ' "> '' 



Shear elastic energy " (2/5)G(7A2 + 12A + 108) G TA^ + 12A + 108 
In terms of the dimensionless variable 



G' 

we have A — 6/(3fc + 1) which, when substituted into Eq. (pi)l . leads to 



fc EE - , (26) 



Bulk elastic energy 2fc 

- n;,2 ( oj, , /I ■ (^'^) 



Shear elastic energy 9^^ + Bfc + 4 

The derivative of the fraction on the right with respect to fc is zero when 2(9fc2 + Bfc + 4) = 2fc(18fc + 8), implying 
h = 2/3. Thus the maximum bulk elastic energy is when k = 2/3. In conclusion, 

Bulk elastic energy Bulk elastic energy ( ^ 2^^ ^ 1 



Shear elastic energy Shear elastic energy 



('-i)^i^^ p"' 



More typically, fc = 5/2 leads to a maximum bulk shear ratio of 20/337 which is roughly 6%. 

One important unsolved problem of condensed-matter physics is the origin of the non-Arrhenius average relaxation 
time of glass- forming liquids, where in most cases one observes an activation energy that increases quite a lot upon 
cooling. One class of theories are the elastic models (recently reviewed in [3| ) according to which the activation energy 
is some linear combination of the instantaneous bulk and shear moduli. As the regards the temperature dependence 
of the activation energy the shear modulus completely dominates (contributing at least 90% [1]). A so-called flow 
event - the jump in configuration space from one to another potential energy minumum - is usually localized in real 
space. Since the surrounding ultraviscous liquid on the short time scale may be regarded as a solid, the molecular 
displacements induced by a flow event may be regarded as those of a point defect (in a disordered solid, albeit). The 
above result sheds light on the nature of far-field flow-event induced displacements by emphasizing "shear dominance" 
- that the shear modulus is much more important than the bulk. 

In an interesting book from 1986 Varotsos and Alexopoulos (VA) discussed point defect energies in solids, as well 
as activation energies for point defect diffusion [19|. VA concluded that these energies can always be written as the 
(isothermal) bulk modulus times some microscopic volume. Thus point defect energies and activation energies always 
scale with the bulk modulus. This result is at variance with the above calculation. In practical terms, of course, 
Poisson's ratio usually does not vary very much, so the bulk and shear moduli are roughly proportional. This means 
that in many cases it is difficult to determine whether the bulk or the shear modulus controls things. Admittedly, 
VA explicitly argue in their book that it is the bulk, and not the shear modulus, which is important, but they do 
not base this on calculations similar to ours (but, e.g., thermodynamic arguments). We do not claim that the above 
calculation disproves VA, but it may encourage other researchers to look into the question. Surely, the idea that many 
important physical properties are controlled by as simple quantities as the elastic moduli deserves to be investigated 
fully. 

This work was supported by the Danish National Research Foundation's Centre for Viscous Liquid Dynamics "Glass 
and Time." 
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